Conductivity and Fano factor in disordered graphene 
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Using the recursive Green's function method, we study the problem of electron transport in a disordered 
single-layer graphene sheet. The conductivity is of order /h and its dependence on the carrier density has a 
scaling form that is controlled solely by the disorder strength and the ratio between the sample size and the cor- 
relation length of the disorder potential. The shot noise Fano factor is shown to be nearly density independent for 
sufficiently strong disorder, with a narrow structure appearing at the neutrality point only for weakly disordered 
samples. Our results are in good agreement with experiments and provide a way for extracting microscopic 
information about the magnitude of extrinsic disorder in graphene. 

PACS numbers: 81.05.Uw, 73.23.-b, 73.50.-h 
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Graphene, a two-dimensional (2D) allotrope of carbon with 
a honeycomb lattice, has attracted a lot of attention due to 
its unusual, Dirac-like, electronic spectrum and its potential 
for an all-carbon based electronics.' While the theoretical 
literature on graphene is already quite extensive)^ the effect 
of disorder on transport properties near the charge neutrality 
point (Dirac point) is still subject to much debate and con- 
troversy. The difficulty in understanding how disorder affects 
Dirac fermions has to do with the fact that electrochemical 
disorder is a relevant perturbation under the renormalization 
group and drives the system away from the weak disorder 
regime.^ Hence, standard perturbation theory fails and one 
has to rely on either non-perturbative methods or numerical 
approaches. Although several recent theoretical works have 
been established that for short-range disorder large graphene 
samples should become insulators^i the current understanding 
of long-range disorder is less clear Some authors^ have ques- 
tioned the existence of a beta function for undoped graphene 
while others^ have proposed a non-monotonic beta function 
for single-valley Dirac fermions and a metal-insulator tran- 
sition. Numerical computations in momentum space^ as well 
as simulations based on a transfer-matrix method*^-^ adapted to 
the single-valley Dirac equation'" found instead a simple scal- 
ing law for the conductivity, a metallic beta function, and no 
indication of a new fixed point. ' ' Furthermore, while the self- 
consistent Born approximation (SCBA) predicts a universal 
(impurity independent) conductivity of 4e^/ (7r/i)r^ percola- 
tion theory finds a smaller value. 

On the experimental side, graphene behaves as a good metal 
with conductivities of the order of e^/h,^* which is inconsis- 
tent with a purely ballistic transport. No sign of strong lo- 
calization has been seen in graphene although the size of the 
samples (a few micrometers) may be smaller than the local- 
ization length in 2D. 

Current experiments indicate that transport is not ballistic, 
but is it really diffusive? Experimentally, the electron mean 
free path iti can be estimated by gating graphene away from 
the Dirac point and using the Drude formula to relate £ti to the 
conductivity a: £tr = tT/[(2e^//i) VrnJ, where n is the car- 
rier density. Typically, £tr is found to be smaller but close to 
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FIG. 1: (Color online) Color maps: [(a)] Long-range disordered po- 
tential in a square graphene sheet (i\f — 161, N — 92, A'o ~ 1, 
7^imp = 0.02, ^ — 4ao); Carrier densities (positive for electron and 
negative for holes) at the Dirac point |(b)[ above it |(c)[ and below 
1(d) I with jj./t — 0, 0.1, and —0.1, respectively. 



the sheet linear size L. Although this estimate has to be taken 
with caution, it supports the notion that cuiTent experiments 
are in the crossover regime between ballistic and diffusive, a 
regime where standard analytical approaches fail. Establish- 
ing a quantitatively understanding of the effect of disorder in 
the transport properties of graphene is of fundamental impor- 
tance for the development of electronic devices where high 
mobilities are required. Therefore, an alternative theoretical 
study is desirable. 

In this paper, we use the recursive Green's function 
method'^ to study numerically the effect of extrinsic long- 
range disorder in the linear electronic transport of graphene 
at low carrier densities. We explore the low temperature 
limit, when intrinsic scattering mechanisms such as electron- 
phonon can be neglected. Our calculations are based on a 
microscopic, tight-binding model of the graphene sheet and 
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performed at and near the Dirac point. After a simple rescal- 
ing, we find that the dependence of the conductivity on the 
carrier density is controlled only by a dimensionless disorder 
amplitude and the ratio of system size to disorder correlation 
length, L/^. This dependence can vary from sublinear to su- 
perlinear, depending upon the disorder strength. Therefore, 
our results provide a way for extracting extrinsic disorder pa- 
rameters from experimental data without making any assump- 
tion about the underlying dynamical regime. We also find that 
disorder has a marked effect on how the shot noise Fano fac- 
tor depends on the carrier density. For samples with square 
aspect ratios and weak disorder, a narrow dip exists at the 
Dirac point, while for sufficiently strong disorder, the Fano 
factor becomes density independent, taking up values larger 
than 1/3 which increase slowly with the disorder 

The nearest-neighbor tight-binding model on a honeycomb 
lattice with site disorder reads H = — X](ij)^N)OI + 
J2i U {vi) \i) where i and label the lattice sites and their 
coordinates, respectively, and t « 2.8 eV is the hopping en- 
ergy. The long-range disordered potential [/(r^) consists of 
A/imp random lattice sites {R,i}. Each of these sites is the 
center of a Gaussian scatterer with a random amplitude C/„ 
taken from a uniform distribution over the interval [—(5,(5]: 

U{ri) = EfiTC^ne-l'-'-i^"l'/2€', where e > ao is the 
range of the potential (ao ~ 2.46 A is the lattice constant). '^ '^ 
A typical realization of the disorder is shown in Fig. The 
concentration of scatterers is nimp = A/imp/A/", where Af de- 
notes the total number of atomic sites. When riimp ^ 1, the 
magnitude of the disorder fluctuations is characterized by the 
dimensionless parameter Kq, which is defined from the impu- 
rity potential correlation function 

where v = \/3actt/2h is the Fermi velocity (notice that 
([/(r^)) — 0). We note that Kq contains information not only 
about the relative magnitude of the potential fluctuations, 6/t, 
but also about the scatterers' range and concentration: A sim- 
ple calculation yields Kq w 40.5 fiiinp(5/0^(C/"o)'*-'^ In the 
continuum limit and using Eq. ([T]i together with the Born ap- 
proximation (BA), one finds that the transport mean free path 
away from the Dirac point is given by if/^ — 2Xf / {nKo), 
where Af is the Fermi wavelength in the graphene sheet 
(Af < ef^^ for flie BA to hold). 

We investigate square graphene sheets with M = l+2L/ao 
andA^ = l+2V3i/ao lattice sites alongperpendicularedges. 
The sheet is connected to ballistic, square lattice source and 
drain leads following standard procedures. '^ '^ For all data 
shown here, the edges parallel to the electron current flow are 
in the metallic armchair configuration'^ (we have found that, 
in the presence of disorder, similar results are obtained for the 
zigzag configuration). Using the recursive method, we com- 
pute the electron retarded Green's function at any lattice site, 
G^{ri;E), where E is the electron energy. The local den- 
sity of states, v{vi\E) = —(l/7r)ImG(ri; S), is obtained and 
can be integrated in energy to yield the local carrier density 
at the Fermi energy, n{vi]^) = /^(j. ) ^{^i'l^)- This is 
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FIG. 2: (Color online) (a) Scaling of the average conductivity with 
sample size for a square graphene sheet at the Dirac point. The 
data corresponds to A/ = 41, 83, 113, 161, 317, 479, ^/ao = 
2, 3, 5, and nimp = 0.02. The lines denotes the best-fits of the form 
cr/(2e^//i) = A\a{L/Q + B io the data, (b) Average conductivity 
versus rescaled sheet lengths, L* = a(A'o) L, with a determined 
from the best-fits found in (a). 



illustrated in Figs. [HElldl where we show the local carrier den- 
sity for a fixed realization of the disorder potential and three 
values of the Fermi energy. The hole and electron puddles 
formed at the maxima and minima of the potential are clearly 
visible. These patterns are reminiscent of the single-electron 
transistor scannings of graphene.^" 

The recursive Green's function method provides a effi- 
cient way to compute the electron's scattering matrix as a 
function of the electron energy, S{E). From 5, one ob- 
tains the transmission t and reflection r matrices, which are 
then used in the computation of the linear conductance Q 
and the shot noise Fano factor F: Q = (2e^//i)Tr(ttt) and 
F = TT{tHr''r)/TT{tH) (the factor of 2 in ^ accounts for 
spin degeneracy). We note that a = Q for a square sheet. 

In Fig. |2l we show the scaling of the conductivity with 
L at the Dirac point for several disorder strengths and cor- 
relation lengths. Our results are in qualitative agreement 
with recent numerical studies using the continuum model^^ 
and are compatible with a metallic beta function of the form 
/5((t) = A {2e^/h) /a, with A w 0.17. In particular, the quan- 
titative agreement with the scaling obtained in Refs. 01^ in- 
dicates that for correlation lengths ^ > 2ao, there is no no- 
ticeable effect of intervalley scattering. Thus, the lattices in 
our study are large enough to reproduce the numerical results 
obtained with the continuous single-valley model. 

In order to explore how the conductivity behaves as a func- 
tion of carrier density n, we compute the latter as a func- 
tion of the Fermi energy fi in the graphene sheet through 
the expression n(/i) = dE ^{E), where i^{E) is the 
sheet global density of states (DOS), which can be readily 
obtained from the energy dependence of the scattering ma- 
trix. By writing the scattering matrix in terms of the re- 
tarded Green's functions,'^ one can show that the Wigner- 
Smith time delay, ty^{E) = —ifi{d/ dE) ln(dct 5), becomes 
Tw(£^) = -2;iTr [ImG"^(ri ;£:)]. We thus arrive at v{E) = 
-{i/TTWL)Tr:[S'^{dS/dE)\. These steps show that v = 
2t-^ / {hW L), as expected from scattering theory. 

In Fig.[3t, we show the average v{E) for different values of 
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FIG. 3: (Color online) (a) Average density of states for three lattice 
sizes (different symbols) and four values of Kq (different colors and 
line types). Notice that L = (A/ - l)ao/2. Only data for £ > 
is shown. The thick solid line represents the case of a clean, infi- 
nite graphene sheet and the other lines are guides to the eyes. The 
number of realizations varies from 100 to 10,000 with ^ = 2ao and 
n-iYnp = 0.02. (b) Shifted density of states for M = 161. The lines 
correspond to best-fits of the form u{E) - iy{0) = Ctal \E/t\°'. 



the disorder parameter A'o and a fixed ^. Notice that the DOS 
does not vanish at the Dirac point, as one expects for an in- 
finite clean graphene sheet. The coupling to the semi-infinite 
leads broadens the electronic energy levels in the graphene 
sheet, an effect that is more pronounced in the vicinity of the 
Dirac point. More interesting however is the effect of disor- 
der. We observe that as Kq increases, the DOS becomes al- 
most flat over a large energy range. In particular, close to the 
Dirac point the DOS grows substantially with Kq, surpassing 
the broadening due to the coupling to the leads. 

Remarkably, we find numerical evidence that the DOS is 
very robust against system size scaling. Figure[3^ shows that 
iy{E) rapidly converges as AI is increased and Kq < 10. An- 
other feature worth noticing is that, for strong enough disor- 
der, i^{E) grows approximately as a power-law with E, as 
shown in Fig. [3j). Our results for the DOS complement the 
analysis of Ref. 3 of the problem of Dirac fermions in the 
presence of a random chemical potential. This kind of disor- 
der was shown to be a relevant perturbation in the renormal- 
ization group sense, possibly taking the system to a new fixed 
point if electron-electron interactions are taken into account.* 

We now turn our attention to transport properties. Our main 
results are collected in Fig. H) As one departs from the Dirac 
point, the counterintuitive feature that the conductivity is en- 
hanced by disorder disappears.— For a fixed disorder, the con- 
ductivity increases with carrier density. By plotting cr as a 
function of the number of electrons contained in a square of 
area we find that the functional dependence is solely con- 
trolled by the disorder strength Kq and the ratio L/(_. Notice 
that cr as a function of rt^^ goes from sublinear (weak dis- 
order) to linear (intermediate disorder) and then to superlin- 
ear (strong disorder) as Kq is varied. All three behaviors are 
seen in the experiments'"* and therefore may be related to the 
strength of the long-range disorder in the samples. For a fixed 
Kq and dilute scatterers, we found little dependence of our 
results on riimp. 

The scaling behavior of cr = <j{n£^'^,Ko,L/£^) works in 
the crossover regime from ballistic to diffusive regimes and 
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FIG. 4: (Color online) Average conductance (left column) and 
Drude transport mean free path (right column) as functions of 
carrier density for several disorder strengths (/("o = 1,2,5), 
lattice sizes (L/ao — 80, 101, 164), and disorder correlation 
lengths i^/ao — 2, 2.5, 3, 3.3, 4, 5). The number of realiza- 
tions varies from 100 to 5,000 depending on Ko and L. The thick 
solid line corresponds to the clean system Fermi wavelength, 
\f ~ 2-y/ n/n. The transport mean free path is calculated using 
ftr = a/[{2eyh)y^. 



departs strongly from the SCBA prediction,'* which is only 
justifiable in the semiclassical diffusive regime, when Aj? <C 



< L. Notice that 



Ai? in Fig. |4]and therefore 



our simulations are far from the semiclassical limit. Although 
the lattices investigated are rather small in comparison to the 
graphene samples probed in experiments, the scaling allows 
us to make a direct comparison with experimental data by us- 
ing realistic values of ^which ranges between 10 and 100 
nm, as shown in Ref. 20) and extrapolating the cr curves to 
large caiTier densities, typically of order 10'^^ cm^^. This in 
turn provides a way for estimating the value of Kq in transport 
measurements dominated by extrinsic disorder. As an exam- 
ple, in Fig. m if we take L/£^ ~ 20 and Kq — 2 (which yields 
an approximately linear cr(n) dependence) and assume ^ = 50 
nm, we obtain a ~ 500 e'^/h, which is just a bit higher than 
the values measured experimentally in high mobility samples. 

In Fig. |5] we present the typical dependence of the average 
Fano factor on the caiTier density and disorder strength for 
square and rectangular sheets. For weakly disordered square 
samples, we observe a narrow dip at the Dirac point followed 
by a saturation of as n is increased. As the disorder in- 
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FIG. 5: (Color online) Average Fano factor as a function of carrier 
density for several disorder strengths and ^ = 2ao. The number of 
realizations ranges from 100 to 5,000. (a) Square sheets (W/L — 
l)withL/^ = 40. (b) Rectangular sheets (W/L = 3.1) with 
L/^ = 26. The dashed lines represent F = 1/3. 



creases, the dip rapidly disappears and F takes a density- 
independent value above 1/3, which is the value expected for 
a normal metal in the diffusive regime.'^ In contrast to the 
behavior seen for the conductivity, we do not observe a clear 
scaling behavior for F, namely, the data do not collapse onto 
distinct universal curves when sorted out by Kq and L/£^ (not 
shown). For short and wide rectangular sheets, the behav- 
ior is quite distinct from the square case. For weak disor- 
der, there is a peak located at the Dirac point, followed by a 
damped oscillation towards a saturation value (the asymmetry 
is due to a band mismatch between the honeycomb and lead 
lattices). These features are compatible with analytical results 
derived for clean junctions with large W/ L aspect ratios.^' As 
the disorder strength increases, the peak widens, the oscilla- 
tions are smeared, and the saturation value of F goes up. For 
large disorder strengths, F becomes again nearly density in- 
dependent. For both geomeries, the behavior we obtain for 
strong disorder appears consistent with recent experiments by 
DiCarlo et al, who found a nearly density independent F for 



several graphene samples.^ However, this behavior departs 
from the results of other numerical simulations, where a satu- 
ration value near 0.3 was found. ^ Notice that the dip (square 
samples) and the peak (rectangular samples) near the Dirac 
point disappear rapidly as the disorder strength is increased 
and therefore can be interpreted as signature of very weak dis- 
order Danneau et al. have recently reported the observation 
of ballistic behavior in their measurements of the Fano fac- 
tor for samples where W/L ^ 1,^'* although the width of the 
peak in their experiment is nearly two orders of magnitude 
larger than the theoretical prediction for clean samples.^' 

In summary, we have performed numerical simulations to 
evaluate the conductivity and the shot noise Fano factor in 
graphene samples as a function carrier density and in the pres- 
ence of long-range disorder We found that the conductivity 
follows universal scaling curves that depend solely on the dis- 
order strength and the ratio of sheet length to disorder coiTe- 
lation length. The Fano factor was found to have a marked 
dependence on the density only for very weak disorder and 
near the Dirac point. For moderate to strong disorder, the 
Fano factor is nearly constant and its value is not universal. 
While our model of long-range disorder applies mainly to sit- 
uations where the strong modulation of the potential in the 
graphene sheet is due to roughness in the substrate, it could 
also be extended to account for the effect of charge impurities 
present in the substrate or on the graphene sheet. The only 
necessary modification is the inclusion of a density-dependent 
correlation length ^ to take into account screening by carriers 
in the graphene sheet. However, this modulation of ^ should 
be weak for the low densities we have been able to study. 
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